A model for coexistent superconductivity and ferromagnetism 
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We explore the various temperature dependencies and thermodynamic quantities of a mean field 
model of a ferromagnetic-superconducting system. The starting point for this model is based on 
an s-wave pairing scheme in the singlet channel of the superconducting state and a spontaneously 
broken symmetry phase in the ferromagnetic state. We show numerically and analytically that a 
state of coexistence reveals itself and is favored energetically over other possible states, and a simple 
phase diagram is developed. Finally, a comparison of the specific heat with experiment is shown. 
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Much work in the past several decades has been done 
studying the interplay between magnetism and super- 
conductivity. Some of the early progress by Abrikosov 

SjEind the work leading to the so-called LOFF model 
yj was directed towards ferromagnetic systems of mag- 
netic impurities coexisting with superconducting order. 
Fay and Appel Q showed how itinerant ferromagnetism 
can promote p-wave superconductivity, and much of the 
theoretical studies was concentrated on developing the 
symmetry of the superconducting order parameter. The 
question of possible s-wave singlet pairing was thought 
possible on the paramagnetic side of the magnetic tran- 
sition, but it was shown even there that ferromametic 
fluctuations destroy it near the Curie temperature [3 • 

More recently, two of the current authors studied a 
weak ferromagnetic Fermi liquid and showed that s-wave 
superconductivity is possible, and that it is the favored 
phase on the ferromagnetic side 0, ■ A possible phase 
diagram was predicted, however, at the time there were 
no materials that demonstrated superconductivity inside 
the ferromagnetic state. Recently a similar considera- 
tion of s-wave superconductivity has been carried out by 
Suhl and Abrikosov in which the ferromagnetism 
is due to localized spins. We instead develop an itinerant 
ferromagnetic model in which the magnetic electrons are 
also the ones responsible for the formation of the Cooper 
pairs. In this work, it must be noted that the BCS state 
is actually the preferred one energetically, as shown in 
ref.^3- However, we are assuming that there already 
exists weak ferromagnetic order, inside which arises su- 
perconductivity. The point is that we are assuming the 
strongest competition to be between the coexistent and 
the ferromagnetic state, with the pure BCS state not 
being a possibility. Moreover, in the microscopic mod- 
els discussed 0, 0, H I3i s-wave superconductivity 
pairing interaction vanishes when the ferromagnetic or- 
der vanishes. With that said, we can consider this an 
effective mean field that has nontrivial results and can 
be shown to coincide with what is seen experimentally. 

The phenomenon of coexistent superconductivity and 
ferromagnetism was then discovered in UGe2 [ul Il2l |. 



and subsequently in ZrZnz [H and URhGe The 
phase diagram of UGe2 is very similar to the one pro- 
posed in refs.j^ Q, however, no superconductivity was 
detected in any of these materials on the paramagnetic 
side, contrary to early predictions. These currrent dis- 
coveries are promoting much more theoretical work to 
explain the superconducting pairing, with most calcula- 
tions based on a p-wave triplet state 0|. Also, it has 
been verified that the magnetic transition in at least 
one of these materials, UGe2, is of first-order below a 
certain temperature. We will not discuss the thermody- 
namics of that transition here, as it will be reserved for 
a future publication. 

In this paper we will present some new and unexpected 
results of an alternative model of s-wave singlet pairing 
developed by Karchev et al.0. In that paper, some 
zero temperature results were calculated, and the finite 
temperature results will be shown here. 

We refer the reader to ref. for the details, and here 
we look at the mean-field Hamiltonian obtained from a 
model Hamiltonian by the standard mean-field proce- 
dure, 
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The diagonalization of this Hamiltonian using a Bogoli- 
ubov transformation yields, 
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The quasiparticle energy dispersion relations are, 
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The final step is to minimize the free energy to produce 
the mean-field equations. This results in a set of two 
coupled equations in M and A that will be solved self- 
consistently below. For M we find, 
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and for A, 
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The two order parameters, M and A, have dependen- 
cies such as M = M{g, J) and A = A(g, J) at T = 0. 
From the numerical solutions of the two coupled equa- 
tions H5I6|I. we show characteristic curves of the order 
parameters, a phase diagram, and the free energy differ- 
ence of different states. We stress that all of the future 
results are derived strictly from the dispersion relations 
and the mean-field equations only, with no other assump- 
tions made about coupling strength limits or small mag- 
netization. 

At T = 0, only one of the fermionic particles con- 
tributes, since E'^ > for all p, but < when p < Pp 
and p > , where 



p| = y2m*fi±m*y^{JM)^ - 4|A|2. 

Therefore, the coupled equations H5I6|I reduce to 
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The only solutions of eqs.® and © which exhibit the co- 
existence of superconductivity and ferromagnetism are in 
the case when JM > 2|A|. However, in theory, there are 
two other solutions: the pure ferromagnetic state with 
A = and the BCS state with M — (and obviously 
the normal, paramagnetic state). We will not consider 
the 'BCS' state here (see introduction). 




FIG. 1: M and A curves at T = as functions of J and g. 
In (a) we see that these order parameters are weak functions 
of J, while in (b) there is strong p-dependence. Note: Here 
and elsewhere, by g we mean 2-K^gN{0), so that these calcula- 
tions are definitely in the weak-coupling regime{gN{0) << 1). 
Also, the magnetization is normalized by n, so that M/n < 1. 



Our concern will be with the case JM > 2|A| as in 
ref.|l7l|. A note is in order here concerning the approxi- 
mations made in refs. |[rRITl| . Analytically, the solutions 
obtained from the coupled equations at T = when JM 
is close to, but larger than 2A are 
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where Aq is the usual energy gap equation in a non- 
magnetic superconducting state: 



Aq — 2ec exp 
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The important thing to notice is that in eq. H10|l the mag- 
netization and the gap are strictly proportional to each 
other in the coexistent state. This is an artifact of the 
approximation used, and we would like to stress that 
the numerical calculations do not make use of this ap- 
proximation. To show this explicitly, we look at Fig^ 
where the coexistent solution with both order parame- 
ters is shown as a function of J and g respectively. 

It is seen from this figure that both order parameters 
go to zero together either at J = 1 or 5 = 0. This is due 
to the self-consistency of the coupled equations. As J 
is increased the order parameters are fairly constant and 
solutions can be found for large values of J. However, we 
are mainly interested in weak ferromagnetism, i.e., when 
the coupling J is close to 1, and we take the value of 
J — 1.001. Using this J we show the order parameters 
as functions of g in (b) of Fig^ This has an interesting 
feature in that the superconducting gap increases from 
g = and then gets suppressed at a higher value of g. 
The magnetization is constant after the gap disappears 
because M is not an explicit function of g, only implic- 
itly through A. This feature of the gap increasing and 
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then being suppressed is non-existent in the approximate 
eqns. . 

In FigI21we have calculated the free energy of the co- 
existent state with relation to the normal state and the 
pure ferromagnetic state (A = 0). In the range where 
the superconducting gap has its largest value, 5 « 2 (see 
Fig^), the free energy of the superconducting ferromag- 
netic state is lower than the two other states, and hence 
is the preferred one. For very small values of g, we see 
a crossover. This behavior will also be evident in the fi- 
nite temperature free energy, where the area around the 
largest value for in the superconducting state yields a 
lower free energy. 

At finite temperatures one can look at the same types 
of dependencies of the order parameters. We solve 
eatns. H5l6|l introducing a debye cutoff in the otherwise 
divergent integrals. The cutoff is small compared to the 
Fermi energy as this is the region where the important 
physics takes place, and varying this cutoff has no quali- 
tative effects on the various quantities calculated below. 
A phase diagram can also be developed which shows how 
the Tc's depend on the couphng constants. Then, ther- 
modynamic properties can be calculated such as the free 
energy and specific heat. 

The temperature dependence of M and A can differ 
widely with respect to changes in the coupling constant 
space. This is why this model is convenient because it is 
adaptable for different cases. By tuning the parameters 
we can see physically realizable results. Figure 13 shows 
some examples. In (a), g = and we see the full mag- 
netization as it approaches its Curie temperature. There 
is no superconducting gap in this case. In the second 
and third graphs, we now turn on the superconducting 
coupling and see that it has dramatic effects on the mag- 
netization. There are two cases. For a small g, the Tc 
for M and A coincide. The superconductivity weakens 
the magnetization and they become proportional and si- 
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FIG. 2: Free energy of various states with respect to the nor- 
mal state(A = M = 0). Notice where the coexistent state has 
its lowest free energy value (see Fig0, i.e., where the gap has 
its largest value. 
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FIG. 3: The temperature dependence of the order parameters 
for different coupling values of g. Here, J — 1.001. The 
temperature, T, is unitless, but can be converted to Kelvin 
with one input parameter, the Fermi energy. For example, 
using this conversion for U Ge2 , the approximate temperature 
is ~ 5r[K]. This gives a Curie temperature of the right order 
of magnitude. 
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FIG. 4: A typical phase diagram for a fixed value of J=1.001 
with the Tc's on the y-axis. Also shown for comparison is the 
Tc dependence of a pure BCS superconductor on g. 



multaneously go to zero. The second case is for a larger 
g where now the Tc's differ. The gap disappears and the 
magnetization falls back onto the full saturation curve, 
as expected from the equations. 

In Fig^a phase diagram is developed from this model. 
It shows the Tc's of each order parameter of the coex- 
istent solution plotted as functions oi g. In (b) is the 
phase diagram of a normal BCS superconductor without 
magnetic ordering and in comparison with (a), it is ob- 
vious that the effect of the magnetization as it gets large 
is to suppress and then kill superconductivity. The su- 
perconducting Tc goes through a maximum and then is 
suppressed eventually to zero. 

As with the zero temperature case, the free energy of 
the coexistent state is lower than that of the pure fer- 
romagnetic state through a broad temperature range for 
certain values of the coupling constants. In Fig[Slwe show 
one case. This particular plot is calculated with g = 1.2, 
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FIG. 5: The y-axis is labelled as the negative of the free energy 
difference of a certain state and the normal state, in arbitrary 
units. T = TcA where the superconducting gap parameter 
vanishes inside the magnetic state, at the point where the 
coexistant— >ferromagnetic transition is. The dashed curve 
does merge smoothly yet abruptly with the solid curve, at 
which point the two lines become indistinguishable. 



hence it is helpful to refer to Fig0Ja) to see which part 
of the phase diagram we are analyzing. In this case, we 
are looking at the zone where the Tc's of the two order 
parameters do not coincide. 

Experiments have indicated that the specific heat dis- 
plays a very small jump (J7Ge2)0| or an undetectable 
one It also shows a very linear temperature 

dependence at low T. This is significant in our model be- 
cause small jumps can be explained by gaplessness, which 
we have with the /3-fermion. Analytically, the specific 
heat jump can be calculated readily to give. 
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FIG. 6: The specific heat jump is much smaller than 
the BCS one, in this case about 15 percent. The 
(rescaled)experimental data is taken from ref. ,19j 



beyond mean field. It is important to realize that the 
question of the nature of the superconducting order pa- 
rameter is still an open one, and we show that a model 
of s-wave pairing is at least qualitatively viable. 

We graciously thank J.R. Englebrecht, Y. Joglekar, 
N.I. Karchev, A.G. Lcbcd, and P.B. Littlewood for use- 
ful discussions. This work was done with the sup- 
port of DOE Grants No. DEFG0297ER45636 and No. 
60202ER63404. 
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The first term on the right gives the usual BCS specific 
heat jump, however the remaining quantity, which is neg- 
ative, gives the decrease from the BCS value. Thus, the 
jump can be tuned depending on the value of the mag- 
netization. In FigEl is the numerical plot of the specific 
heat jump at the superconducting transition inside the 
coexistent state. The specific heat is calculated with a 
value of 5 = 1.1 which corresponds to 'strong' supercon- 
ductivity, i.e., where Tca is large. This can be seen in 
the phase diagram, Fig^] Also shown for comparison are 
some experimental data points for UGe2 and it is clear 
that the numerical calculations model well the qualita- 
tive behavior of the specific heat. In that experiment, 
ref.|0|, the specific heat is measured at 1.13 GPa where 
superconductivity is maximum. 

In this paper, a mean field solution of s-wave pairing 
was investigated and shown to give nontrivial results and 
can be the starting point for more detailed calculations 
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